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1 Introduction 



The quark model practically describes all phenomena in hadron physics very successfully. 
This model predicts the existence of hadrons containing two heavy quarks and one can 
estimate their masses in this framework [1]. The first experimental observation of the doubly 
heavy baryons was reported by SELEX Collaboration [2], where the double charmmed 
and baryons were found. This observation was confirmed by measurement of a different 
weak decay mode [3,4]. The experimental discovery of the doubly heavy baryons and study 
of their properties constitutes one of the main directions of the physics program at LHC. 
(for a review on doubly heavy baryons, see for example [5]). 

In this prospect, it would be interesting to present reliable theoretical predictions on 
properties of these baryons. The masses of the doubly heavy baryons have been estimated 
in the frameworks of the quark [6-8] and MIT bag [9] models. In order to study the 
properties of these hadrons in a model independent way, the QCD sum rules method [10] 
based on QCD Lagrangian is one of the most reliable approaches. In the present work, we 
calculate the masses and residues of the doubly heavy spin-1/2 baryons within QCD sum 
rules formalism. We consider both symmetric and anti-symmetric currents with respect to 
the exchange of heavy quarks defining the baryons under consideration in their most general 
forms. It should be noted that the masses of doubly heavy baryons with spin-1/2 have been 
calculated within the same framework using the Ioffe current in [11, 12] and general current 
in [13]. However, the obtained expressions for mass sum rules in [13] at Ioffe current case 
does not exactly reduce to those presented in [11, 12] for Ioffe current. The masses of the 
doubly heavy baryons with spin-3/2 have also been investigated in [13] and [14]. 

The plan of this work is as follows. In next section, the QCD sum rules for the masses 
and residues of the doubly spin-1/2 baryons are obtained. Section 3 is devoted to the 
numerical analysis of the sum rules for physical quantities under consideration. This section 
encompasses also our comparison of the obtained results with those existing in the literature 
as well as our concluding remarks. 

2 Sum rules for the masses and residues of the doubly 
heavy baryons with spin-1/2 

Before presenting the detailed calculation of sum rules for the masses and residues of the 
doubly heavy baryons let us discuss the ground states of these baryons in quark model. 



1 



For the ground state of doubly heavy baryons with the identical heavy quarks, Eqq and 
figg, the pair of heavy quarks form a diquark with total spin of 1. Adding then the spin- 
1/2 of light quark, we have two states with total spin-1/2 or spin-3/2. The baryons with 
(without) star stand for spin 3/2 (1/2) doubly heavy baryons during the text. For these 
states, the interpolating current should be symmetric with respect to the exchange of the 
heavy quarks fields. For the states containing two different heavy quarks, in addition to the 
previous case, i.e., total spin of diquark equal to one, the diquark can also have the total 
spin zero, which leads to the total spin 1/2 of these states. Obviously, the interpolating 
current of these states (usually these states are denoted by E' bc and fl' bc ) are anti-symmetric 
with respect to two heavy quarks fields. In the present work, we deal only with spin-1/2 
doubly heavy baryons. 

After these preliminary remarks, let us come back to our main problem, i.e., calculation 
of the masses and residues of the doubly heavy baryons. The interpolating currents here 
play the role of the wave functions in the quark model. The general expressions of the 
interpolating currents for the spin-1/2 doubly heavy baryons in their symmetric and anti- 
symmetric forms can be written as: 

V s = -^e abc ^(Q aT Cq b ) l5 Q' c + {Q' aT Cq b ) lb Q c + f3(Q aT C l5 q b )Q' c + l3{Q' aT C^q b )Q c 
V A = -^=e abc ^2(Q aT CQ' b ) l5 q c + (Q aT Cq b ) l5 Q' c - (Q' aT Cq b ) l5 Q c + 2^Q aT C l5 Q' b )q c 

+ m aT cid)Q ,c - m aT cid)Q c ^ (2.1) 

where (3 is an arbitrary auxiliary parameter. The case, (3 = — 1 corresponds to the Ioffe 
current. Here C stands for the charge conjugation operator, T denotes the transposition, 
a, b, and c are the color indices; and Q and q correspond to the heavy and light quarks 
fields, respectively. The interpolating current with the light quark u or d corresponds to 
the ^QQ q , and with s to the Qqq q baryons, respectively. Here we would like to note once 
more that in the symmetric part, both heavy quarks may be identical or different, but in 
the anti-symmetric part two heavy quarks must be different. 

For calculation of the masses and residues of the doubly heavy baryons within the QCD 
sum rules, we consider the following correlation function: 

U s ^ A \q)=t J rf 4 xe l9:c (0|r{^ 5(A) (^)^ (A) (0)}|0), (2.2) 
where q is four-momentum of the doubly heavy baryon. This correlation function can be 



2 



written in terms of two independent structures 

n^)(g) = ^(g 2 ) + f /n^(g 2 ), (2.3) 

where U stands for unit matrix. 

The above mentioned correlation function can be calculated in two different manners. 
From one side, it is calculated in terms of hadronic parameters called physical or phe- 
nomeno logical part. On the other side, it is evaluated in terms of quarks and gluons and 
their interactions with each other and QCD vacuum called theoretical or QCD side. Equat- 
ing these two different representations of the correlation function to each other under the 
flag of the quark-hadron duality assumption, we get sum rules for the physical quantities 
under consideration. To suppress the contributions coming from the higher states and con- 
tinuum and obtain contribution of the ground state, we apply the Borel transformation and 
continuum subtraction after performing the Fourier integrals. These procedures bring two 
auxiliary parameters so called Borel mass parameter and continuum threshold, which we 
shall also find their working region to numerically analyze the obtained sum rules in the 
next section. 

Saturating the correlation function by a complete set of hadronic states with the same 
quantum numbers as the interpolating current and isolating the ground state baryons in 
phenomenological side, we have 

u s(A ){) = (0\v s ^(0)\B( q ))(B( q )\f^(0)\0) + {2A) 

q 2 -m 2 B '"' 

where dots refers to the contribution of the higher states and continuum. The matrix 
element {0\rj(0)\B(q}) for spin-1/2 baryons is determined as 

(0\ V (0)\B(q,s)) = X B u(q,s), (2.5) 

where \ B is the residue and u(q, s) is the Dirac spinor with s stands for the spin. Putting 
the above equations all together, and performing summation over spins, we get the following 
final representation for the phenomenological side: 



U S(A)^ = x2 B s ( A)(<i + m B s ( A)) + 



2 (2.6) 

Q m B S(A) 



where we have only two independent Lorentz structures ^ and U discussed above. 

In QCD side, the correlation function is calculated in deep Euclidean region with the 
help of operator product expansion (OPE). Using the Wick theorem and performing all 
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contractions of the quarks fields for the symmetric part, we get the following expression in 
terms of heavy and light quarks propagators: 

U s (q) = iAe ahc e a , Vc , j rf 4 xe^(0 | { - -y B Sg S?< S# n - 1,S$ S^' S<$' lb 



+ lb S% lb Tr 



Qab' Qlba' 



cab' Qlba' 
D Q' 



+ p(- i 5 s$' l5 s' q ba 's a Q i - l5 s$ l5 s' ba 's a Q c ' - s$'s' q ba ' l5 s a Q c ; l5 

- S$S' ba ' l5 S a Q c ' l5 + l5 S#Tr [sgjtSf] + S^l^Tr [s ab> S'^' lb ] 



+ l5 S c Q c Tr S$nS, 



-tfba' 



+ 



S% ' lb Tr[s ab ; 'S' ba ' 75 ]) 
/3 2 ( - SgnSfnStf - S$ l5 S lba ' l5 S a Q c ' 



+ S c n c ,Tr 



75^Q 75 



S c Q c Tr 



^'75^75])} I 0) 



(2.7) 



where S' = CS T C. In the case of Q 7^ Q', the constant A in the above equation takes the 
value A — |, while when Q = Q' we have A = 1 as a result of extra contractions between 
the same quark fields. For the anti-symmetric part we have 

H A (q) = l -e abc e a , b/& fd 4 xe^(0\ fasg S%?' S bc ' l5 + l5 S$ S' ba ' S a Q i l5 



6 

2 l5 SfS'fS b Q c ', l5 + 2 7s Sf 'S' bb 'S a Q c '75 + 4 75 ^ c ' 75 Tr 



oaf/ Qlba' 



oaf/ o/6a' 



/af/ c^c' 
1 



+ 7 5^'75Tr[^'^'j + 7 5^ c '75Tr 

+ (275^75^' + 75^'75^°'^' - 275 ^75^^, 

+ lsS<*: ib S> q ba 'S%' - 2 l5 Sf l5 S% b 'S b Q c ', + 2 l5 Sf l5 S' bb 'S a Q c ' 

+ 2S$S^' l5 S bc ' l5 + SgSfTtSgjt - 2S#Sg 1 s&f 1 s 



+ SgsF-nSgTt - 2SfS% b ' l5 S b Q c ; i5 + 2S?S'$ 7 5^ c 7s 



+ 4 75 ^ cc 'Tr Sj 76 5g 



+ ^7sTr 



rm£>' Qlba' _ 
^ 75 



+ 4,S g cc 75 Tr 
+ 75^q 'Tr 



oaf/ Qlba' 
°Q °Q' 75 



Qab', Qlba' 

b Q 75^g 



c<ab „ Qlba 
b Q' 75>Jg 



+ -Y S S$Tr 
S%' l5 Tr[s$S' ba ' l5 ]) 



+ 



P 2 (2S§' l5 S%?' l5 S bc ' + Sg-YtSF-YsSg - 2S# l5 S'f l5 S bc ' 



+ S$<kS™<kS% - 25^ 75^"° 75^' + 2S?tsS%tsS% 



-iiba 



lab' 



+ 45 g cc Tr 



nba' , QiaO', 
^O'I^q 75 



/af/ 



^ a 75^75 



/af/ 



^ a '75^'75])} I 0). 



(2.8) 



For calculation of the QCD side, we need to know the explicit expressions of the light 
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and heavy quarks propagators. Their expressions in coordinate space are presented in the 
Appendix. 

The coefficient of any structure in Eq. (2.3) in QCD side can be written as the following 
dispersion integral: 

ro s(A) (s) 



n^V) = J '-f^ds, (2.9) 

where pf^ are spectral densities which can be obtained from imaginary parts of the cor- 
relation functions, i.e., 

pf^( S ) = i/m|nf^( S )|, (2.10) 

where subindex 1(2) corresponds to the coefficient of the structure i{U). 

Our main task now is to calculate these spectral densities. For this aim, we write integral 
representation of the Bessel functions in Euclidean space, 



K v (m Q y/=tf) If dt ^ (f+ fj) ( - 

where x\ = —x 2 . We also use the Schwinger representation for the terms containing 
to write them as exponential form after Wick rotation, i.e., 

—— = —— / d\\ n - 1 e-* x B. (2.12) 

Then, we perform the Gaussian integrals d D XE and make dimensional regularization. After 
lengthy calculations for the spectral densities, we get the results presented in the Appendix. 

As has already been noted, the mass sum rules for the heavy baryons are also calculated 
in [11, 12] and [13]. Here we should remind that although the general form of the inter- 
polating current is used in [13] similar to our case, the results are presented only for the 
symmetric current case and the anti-symmetric current with respect to the heavy quarks is 
not analyzed. On the other hand, the results presented in [11] and [12] are obtained only 
for the Ioffe current (ft = —1). Here we would like to compare our results that are obtained 
using the symmetric current with those presented in [13] as well as our predictions in the 
case of special Ioffe current with those of for instance [11] in details. 

Comparing our results on the spectral densities for the symmetric current case with 
those given in [13], we see that as far as the $ structure is considered, the perturbative 
parts coincide, but there is a sign difference in the terms containing (qq). In the case of 
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the structure U, the situation is reversed, i.e., the signs of the perturbative parts in two 
works are different while the terms containing (qq) have the same sign. Moreover, the terms 
containing the quark condensates multiplied by m q factor also show the same differences. 
For both structures, the coefficients of the operator ml(qq) in our case are totally different 
from those presented in [13]. 

Comparing our results for the Ioffe case and symmetric current with the results of [11], 
we observe that the coefficient of the mQiriQ' factor in spectral density pf(s) in [11] contains 
an extra multiplying factor of 4. Moreover, the quark condensate term is multiplied by the 
factor 5/8, while in our case this factor is 1/2. 

Now let us compare the results for the spectral density pf(s). In our case, there is no 
any term multiplied by only the light quark mass m q (without condensates), while in [11] 
there exists such a term (first term in Eq. (26) in [11]). When we check the term containing 
the quark condensate multiplied by rnQfUQ/, our result is half of that given in [11]. 

In the anti-symmetric current case, we obtain the following differences among our results 
compared to those given in [11]. 

For the spectral density Pi(s): 

• the coefficient of the term mQUiQi in [11] contains an extra factor of 2. 

• The quark condensate term in our case is multiplied by 1/4, but it is multiplied by 
-5/32 in [11]. 



• We have the term, 



— {qq) (1 - a)m Q + am^j . 



12 

in our result, which is totally absent in [11]. 



• Also, the term 



1 Tfl 

T7T n(,P m Q' + Ctm Q)-l 

lb ap 



is absent in [11]. 
When the spectral density p%(s) is analyzed, we see that 
• the term, 

appears in our case, but is absent in [11]. 
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• The coefficients of the operators (qq) and m^(qq) in [11] are also different compared 
to our results. 

Our final task in this section is to equate the coefficients of the selected structures from 
both physical and QCD sides to obtain QCD sum rules for the mass and residue of the 
doubly heavy baryons under consideration. After performing Borel transformation with 
respect to the variable-g 2 as well as continuum subtraction to suppress contribution of the 
higher states and continuum; and using quark-hadron duality assumption we get 

„ - m2 B S(A) f S _ s 

A s s(A)e = / dsp^ J (s)e^. 

J(n 



„2 



l{m Q +m Q ,) 2 
rso 



^ 2 B s(A)m B s(A)e & = / dsp2 ( - A \s)e^ 7 , (2.13) 

J(m Q +m Q ,) 2 

where, M and sq are Borel mass parameter and continuum threshold, respectively. Obvi- 
ously, the mass sum rule for the baryons under consideration is obtained eliminating the 
residue in each of the above equations. This is possible applying derivative with respect to 
the — -pj to both sides of each equation in Eq. (2.13) and dividing by themselves. As a 
result, we get 

= ~P ds^(s)^ ' <2 ' 14) 

where, % can be either 1 or 2. 



3 Numerical results 

In this section, we present our numerical results on the mass and residues of the doubly 
heavy spin-1/2 baryons. For the heavy quark masses, we use their MS masses, which are 
given as, m c (m c ) = (1.28 ±0.03) GeV, m b (m b ) = (4.16 ±0.03) GeV [15] and m s (2 GeV) = 
(102 ± 8) GeV [16]. The values of the quark condensates are taken as uu(l GeV) = 
dd(l GeV) = -(246^ MeV f [17], (ss) = 0.8(uu> and m 2 = (0.8 ± 0.2) GeV 2 . It should 
be noted here that in our calculations we neglect the gluon condensate contributions, since 
their influence to the total result is quite small due to the fact that we have two loops in 
perturbative part and radiation of gluons from heavy quark propagators contains an extra 
a s and I/tuq suppression. This expectation is confirmed by explicit calculation in [12], 
and it is shown that the gluon contribution constitutes only about 2% of the total result. 
Therefore, gluon condensate contribution can safely be neglected. Here we should also 
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stress that the two different structures under consideration give approximately the same 
results, hence in the following, we present the results extracted for the structure </[. 

The sum rules for physical quantities contain also three auxiliary parameters, namely 
Borel mass parameter M 2 continuum threshold s and general parameter j3 entered the 
general spin-1/2 currents. We shall find their working region such that the masses and 
residues are practically independent of these parameters according to the standard criteria 
in QCD sum rules. 

The continuum threshold s is not totally arbitrary but depends on the energy of the 
first excited state. As we have not enough information about the first excited states of 
these baryons, we take the value of the continuum threshold to be in the intervals sq = 



(108 - 120) GeV 2 for bb, s = (45 - 56) GeV 2 for be and s = (13 - 21) GeV 2 for cc baryons 



in accordance with the baryons containing a single heavy quark. Note that there is an 
approach which introduces the effective thresholds, i.e., it is assumed that the continuum 
threshold sq is dependent on Q 2 [18]. But in the present work we will follow the standard 
procedure, i.e., So is independent of Q 2 . Our numerical results show that in the presented 
intervals for different doubly heavy baryons, the physical quantities under consideration 
weakly depend on the continuum threshold. 

The upper bound on the Borel mass parameter M 2 is obtained from the condition that 
the pole contribution is larger compared to continuum and higher states. For this aim, we 
consider the ratio, 



which describes the relative contributions of the continuum and pole. Demanding that 
R > 1/2 (i.e., the pole contribution exceeds the contributions of the higher states and 
continuum), we obtain the maximum values of M 2 for different channels as follows: 



The lower bound on M 2 is determined from the condition that the perturbative contri- 
bution should be larger compared to the nonperturbative contribution. As a result of the 



R 




(3.15) 




6 GeV 2 
9 GeV 2 
15 GeV 2 



(at y/s^ = 4.2 GeV), for S cc and VL CC 
(at y/so = 7.5 GeV), for E bc and Q bc 
(at = 10.9 GeV), for E bb and Vl bb . 



(3.16) 



8 




M 2 (GeV 2 ) 



Figure 1: The dependence of the mass of the E bb baryon on M 2 at the fixed value of 
^ l /s^ =10.9 GeV and at six different values of the parameter (5. 

analysis of this restriction we get, 

{3 GeV 2 , for E cc and VL CC 
6GW 2 ,forS 6c andft 6c (3.17) 
10 GeV 2 , for E bb and Q bb , 

at aforementioned values of ^/s^. 

Our numerical analysis shows that in these "working regions" of M 2 , the perturba- 
tive contributions also exceed the nonperturbative contributions. For example, for the E bb 
baryon, at = 10.9 GeV and at M 2 = 11 GeV 2 the contribution form perturbative 
part constitutes about 70% of the total result, while the higher states and nonperturbative 
contributions constitute about 30% of the total result. Similar results are observed for all 
considered baryons. 

As an example in Fig. (1) we present the dependence of the mass of the E bb baryon 
on M 2 at the fixed value of y/s^ = 10.9 GeV and at six different values of the parameter 
(3. From this figure, we see that the mass of the E bb baryon exhibits a good stability with 
respect to the variation in M 2 . 

Having determined the working region of M 2 and the values of ^/so, our final attempt 
is to find the working region of the arbitrary parameter (3. Depicted in Figure (2) is 
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Figure 2: The dependence of the mass of the E bb baryon on (3 at the fixed value of = 
10.9 GeV and at six different values of M 2 . 

the dependence of the mass of the baryon on /3. We see from this figure that when 
(3 > 2 and /3 < —2 the mass of the S^, baryon depends very weakly on j3 and we obtain 
m Sbb = 9.96 ± 0.90 GeV at M 2 = 11 GeV 2 and ^ = 10.9 GeV. 

Performing similar analysis for other baryons by using the above-mentioned regions for 
the auxiliary parameters, we obtain the numerical values of their masses which are presented 
in Table I. 
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Baryon 


M 2 


v 7 ^ 


This work 


[11] 


[12] 


[13] 


[6] 


[19] 


Exp [20] 




11.0 


10.9 


9.96(0.90) 


9.78(0.07) 


10.17(0.14) 


9.94(0.91) 


10.202 


— 


— 


&bb 


11.0 


10.9 


9.97(0.90) 


9.85(0.07) 


10.32(0.14) 


9.99(0.91) 


10.359 


— 


— 


^bc 


8.0 


7.5 


6.72(0.20) 


6.75(0.05) 


— 


6.86 


6.933 


7.053 


— 


&bc 


8.0 


7.5 


6.75(0.30) 


7.02(0.08) 


— 


6.864 


7.088 


7.148 


— 


' — 'CC 


5.0 


4.6 


3.72(0.20) 


4.26(0.19) 


3.57(0.14) 


3.52(0.06) 


3.620 


3.676 


3.5189(0.0009) 


^cc 


5.0 


4.6 


3.73(0.20) 


4.25(0.20) 


3.71(0.14) 


3.53(0.06) 


3.778 


3.787 






8.0 


7.5 


6.79(0.20) 


6.95(0.08) 






6.963 


7.062 




tt'bc 


8.0 


7.5 


6.80(0.30) 


7.02(0.08) 






7.116 


7.151 





Table 1: The mass of the doubly heavy spin-1/2 baryons (in units of GeV) at j3 — ±2. 



For comparison, we also present predictions of some other theoretical papers as well 
as existing experimental data on mass of S cc doubly charmed baryon. The errors in the 
values of the present work also belong to the uncertainties in determination of the working 
regions for different auxiliary parameters, as well as uncertainties in the values of the input 
parameters. From Table 1, We observe that different approaches predict very close results 
on the masses of the doubly heavy baryons. Note that although there are some discrepancies 
between our results on the expressions of the spectral densities with those of the [11, 12] and 
[13] discussed in the previous section, considering the uncertainties of the values presented 
in this Table, there is a good consistency between our numerical values and those existing 
in [11,12] and [13] which apply the same approach to calculate the masses of the doubly 
heavy baryons. 

We also present the numerical values for the residues of the doubly heavy baryons in the 
present work together with the existing prediction of [11] and [12] in Table 2. The values 
presented inside the parenthesis show the uncertainties of the results on the residues. In 
obtaining the values of residues of the doubly heavy baryons, we use the values of M 2 and 
y/so that are presented in Table 1, at (3 — ±2. It follows from Table 2 that our predictions 
on the residues are quite close to the results of [12], except the values of residues for 
and Qbb- But our predictions and the predictions of [12] on the residues, are considerably 
different compared to the ones given in [11]. 



Baryon 


Present work 


[11] 


[12] 


^bb 


0.44(0.08) 


0.067 -J- 0.057 


0.252(0.064) 


&bb 


0.45(0.08) 




0.311(0.077) 


2 be 


0.28(0.05) 


0.046 ^0.021 




^bc 


0.29(0.05) 






^cc 


0.16(0.03) 


0.042 + 0.026 


0.115(0.027) 


^cc 


0.18(0.04) 




0.138(0.030) 


"6c 


0.30(0.05) 






V'bc 


0.31(0.06) 







Table 2: The residues of the doubly heavy spin-1/2 baryons (in units of GeV 3 ). 
At the end of this section we would like to mention that according to Eq. (2.13), there 
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is another way for determination of the masses of the heavy baryon, i.e., 

mw "W^P' (3 - 18) 

Our numerical calculations show that, using this formula, the central values of the masses 
presented in Table 1 are changed maximally 5%. 

In summary, we have calculated the masses and residues of the doubly heavy spin- 
1/2 baryons considering both symmetric and anti-symmetric currents with respect to the 
exchange of heavy quarks in their most general forms in the framework of QCD sum rules. 
Our results are in a good consistency with the predictions of the other theoretical papers 
as well as the experimental data on the mass of the doubly-charmed S cc baryon. We hope 
that the LHC will provide possibility to study the doubly heavy baryons in near future. 
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4 Appendix 



In this appendix, we present the expressions of the light and heavy quarks propagators as 
well as the spectral densities. 

In our calculations, we have used the following expression for the light quark propagator: 



The heavy quark propagator in an external field is given as 



Sq(x) 



m 



4tt 2 



d A k 



x/ 



— i 
i 

a lkx I du 
o 



UmQ G^{ux)a, u + -^- V2 x,G^ lu 



2(m 2 - k 



2\2 



1717 



(4.20) 



where K 1 and K 2 are the modified Bessel function of the second kind. 

The explicit expressions for the spectral densities obtained after procedure mentioned 
in the body text are given as: 



A 



/ / dtpdr/l 3/x 



1287T 4 

6(-l + /3 2 )m q (rjm Q + ipm Q < 



5 + (3(2 + 5/3) n + 2(-l + V + ?7)(-l + /3) 2 m Q m Q 

-l + 0)^[5 + /3(2 + 5/3) 



+ Mf d^- 



16tt 2 



V>r> 



m 



+ 3(-l + /3 2 ) (-1 + ii)m Q - ijm Q/ 



(4.21) 



+ 



3-077 



5 + /3(2 + 5/3) /i + 2(-l + /3) (-1 + ^ + ^(13 



+ ll/3)m Q m Q/ - (1 + 5/3)m q (r]m Q + ipm Q < 



(gg) 
96tt 2 



AH -3(-l + ^[5 + /3(2 + 5/3)]m g 



+ (-1 + /?)(! + 5/3) [(-l + ^)m Q -^m Q ,]|, 



(4.22) 



14 



p!(s) = 128^7 / #^3^ 

" tymin ^ '/mm I 



3^(-l + (3 2 )fim Q , + m Q 3r){-l + (3 2 )n 



- 2 



5 + /3(2 + 5/3) 



m q m Q , 



+ 



A(qq) 
32tt 2 



dM - 



3m n 



+ V - 2s + 2 5 + (5(2 + 5/3) m Q m Q , + 6(-l + f3 2 )m q (-1 + ip)m Q - ipm Q 



(4.23) 



where, 



1 



256tt 4 



tpmax Pf]max 



dipdr)< ji 



6[5 + (3{2 + 5f3)]m q m Q , 



if)(-l + /3)(1 + 5/3)/im Q , + m Q 7j(-l + /3)(1 + 5/3)/i 

/qq\ Fornax 



1927T 2 



;-i + ^(-i + /3)(i3 + n/3) 



3mp 



+ 4/x' - 2s + 6 5 + 0(2 + 5/3) m Q m Q / + 2(-l + /?)(! + 5/3)m g (-1 + ^)m Q - ipm Q . 



+ 2 m o(l-/?) 2 ^ (4-24) 



// 



/i = 



if) + T) 

rrin m) 



- s, 



T]max 
Ipmax 



if) 1 — if) 

2 

Q' 



ipm 2 , 



sif) — m 2 ' 



1-ij), 
1 



2s 
1 

2^ 



s + itLq — itLq, — \/(s + m 2 n — m 2 ni ) 2 — Am 2 n s 



s + m 2 Q -m 2 Q , + y / (s + m|-m 2 3 ,) 2 -4m 2 3 sj . (4.25) 
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